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Os ■ ABSTRACT 

o 
o 

O ■ We present the first multi-epoch spectropolarimetry of a type II plateau supernova 

^ ! (SN II-P), with optical observations of SN 1999em on days 7, 40, 49, 159, and 163 



after discovery. These data are used to probe the geometry of the electron-scattering 



O \ atmosphere before, during, and after the plateau phase, which ended roughly 90 days 

'^ \ after discovery. Weak continuum polarization with an unchanging polarization angle 

f^ ; [e ^ 160°) is detected at all epochs, with p ^ 0.2% on day 7,pi^ 0.3% on days 40 and 49, 

and p ~ 0.5% in the final observations. Distinct polarization modulations across strong 

^ ' line features are present on days 40, 49, 159, and 163. Uncorrected for interstellar 

j^ . polarization (which is believed to be quite small), polarization peaks are associated 

with strong P-Cygni absorption troughs and nearly complete depolarization is seen 

across the Ha emission profile. The temporal evolution of the continuum polarization 

and sharp changes across lines indicate polarization intrinsic to SN 1999em. When 

modeled in terms of the oblate, electron-scattering atmospheres of Hofiich, the observed 

polarization suggests an asphericity of at least 7% during the period studied. The 

temporal polarization increase may indicate greater asphericity deeper into the ejecta. 

We discuss the implications of asphericity on the use of type II-P supernovae as primary 

extragalactic distance indicators through the Expanding Photosphere Method (EPM). 

If asphericity produces directionally dependant flux and peculiar galaxy motions are 

characterized by cr^^^^ = 300 km s~^, it is shown that the agreement between previous 

EPM measurements of SNe II and distances to the host galaxies predicted by a linear 

Hubble law restrict mean SN II asphericity to values less than 30% (3(t) during the 

photospheric phase. 



Subject headings: distance scale — polarization — supernovae: individual (SN 1999em) 
— techniques: polarimetric 



1. INTRODUCTION 

The assumption of spherical symmetry in core-collapse supernova (SN) explosions is under 
assault. On the observational front, high- velocity "bullets" of matter in SN remnants (e.g., Taylor, 
Manchester, & Lyne 1993), the Galactic distribution (Lyne 1998) and high velocities of pulsars 
(Strom et al. 1995; Cordes & Chernoff 1998), the aspherical morphology of many young SN rem- 
nants (Manchester 1987; see, however, Gaensler 1998), and the asymmetric distribution of material 
inferred from direct speckle imaging of young supernovae (SNe; e.g., SN 1987A, Papaliolios et 
al. 1989; see, however, Hoflich 1990) collectively argue for asymmetry in the explosion mecha- 
nism and/or distribution of SN ejecta. Moreover, recent advances in the understanding of the 
hydrodynamics and distribution of material in the preexplosion core (Bazan &: Arnett 1994; Lai 
& Goldreich 2000) coupled with results obtained through multidimensional numerical explosion 
models (Burrows, Hayes, & Fryxell 1995) imply that asphericity may be a generic feature of the 
explosion process (Burrows 2000). It has even recently been proposed that some core-collapse SNe 
produce gamma-ray bursts (e.g., Bloom et al. 1999) through the action of a jet of material aimed 
fortuitously at the observer, the result of a "bipolar" , jet-induced, SN explosion (Khokhlov et al. 
1999; MacFadyan &: Woosley 1999; Wheeler et al. 2000). Such ideas shake the theoretical founda- 
tion of the conventional spherical neutrino-driven SN explosion (Colgate & White 1966), and beg 
for a direct observational probe of early-time SN geometry. 

Shapiro &: Sutherland (1982) first proposed that polarimetry could provide a powerful tool with 
which to study of young SN geometry. The most basic question is whether the SN is round. Because 
electron scattering dominates the opacity of SNe during the photospheric phase (Branch 1980), any 
continuum radiation that emerges at an angle to the surface normal is linearly polarized. For a 
spherical source that is unresolved (true for all young, recent SNe) the electric vectors cancel exactly, 
resulting in zero net polarization. Aspherical sources, however, will yield residual polarization 
due to incomplete cancellation; see Figure 1 of Leonard, Filippenko, & Matheson (2000b). In 
general, linear polarizations of ~ 1% are expected for moderate {^ 20%) SN asphericity. The exact 
polarization amount varies with the degree of asphericity, as well as with the viewing angle and the 
extension and density profile of the electron-scattering atmosphere (Hoflich 1991a [H91]; Hoflich 
1995). 

SNe come in many varieties, with the main classification based on the presence (type II) or 
absence (type I) of hydrogen in the spectrum near maximum light (see Filippenko 1997 for a recent 
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review). Except for SNe la, whose progenitors are carbon-oxygen white dwarfs that accrete matter 
from a companion and undergo thermonuclear runaway, all other SN types (II-P, II-L, Iln, lib, lb, 
Ic) are believed to result from the explosion of evolved, massive stars (initial mass > 8 — 10 Mq) 
that have reached the ends of their nuclear-burning lives. These core-collapse events present a 
smorgasbord of spectral and photometric properties, leading to a proliferation of subclassifications. 
A consensus is emerging, however, that much of this diversity is due to the state of the progenitor's 
hydrogen envelope (helium as well in some cases) and immediate circumstellar environment (CSE) 
at the time of the explosion. 

"Normal" core-collapse events are thought to result from isolated stars with thick hydrogen 
envelopes (generally of several solar masses) intact at the time of the explosion. Their light curves 
show a distinct plateau (hence the moniker "SN II-P"), resulting from an enduring period (some- 
times as long as 150 days) of nearly constant luminosity as the hydrogen recombination wave recedes 
through the envelope and slowly releases the trapped, shock-deposited energy. The spectral and 
photometric evolution of SNe II-P has been successfully reproduced by computer simulations of 
exploding red supergiant stars (e.g., Eastman, Schmidt, & Kirshner 1996; Weaver & Woosley 1993; 
Imshennik & Nadyozhin 1974; Eastman et al. 1994; Popov 1993). The basic model includes a source 
of thermal photons that electron-scatter their way out of the envelope (Branch 1980), with line fea- 
tures imprinted onto the continuum by resonance scattering in the envelope (see, e.g.. Branch et al. 
2000). The line features in photospheric-phase SNe are characterized by "P-Cygni" profiles, with 
blueshifted absorption troughs with emission centered near the rest wavelength. These character- 
istics are common for lines formed in expanding atmospheres. The general situation is illustrated 
by Figure 1. The absorption component results from continuum photons being scattered out of the 
line-of-sight (1-o-s) by gas approaching the observer in region 1. Photons scattered into the 1-o-s by 
gas in regions 2 and 3 produce the emission component centered near the line's rest wavelength, 
Ao. 




Observer 



Fig. 1. — Schematic representation of an expanding SN atmosphere during the photospheric phase. 
Specific regions referred to in the text are indicated. 



The observed properties of SNe II-P are consistent with those expected for a massive star 
exploding into the near- vacuum of the ambient interstellar medium. Some massive stars, however, 
cast off much of their hydrogen envelope prior to exploding, perhaps through stellar winds or due 
to the influence of a companion. It is the explosions of these stars that are believed to produce 
SNe II-L ("L" for their linearly declining light curves, lacking a plateau) or SNe Iln ("n" for the 
presence of relatively narrow emission lines lacking the customary P-Cygni absorption components 
[Schlegel 1990]). Much of the spectral peculiarity of these events likely results from the interaction 
between the SN and the dense CSE (Chugai 1997; Branch et al. 2000) in which the star explodes. 
Progenitors that have lost their entire hydrogen envelope well before the explosion produce SNe lb 
(Type I due to the absence of hydrogen in their spectra), and those that have been stripped of 
both hydrogen and helium produce SNe Ic. Stars that explode with just a low-mass outer layer of 
hydrogen remaining are placed into the transitional category "SN lib" since the spectrum shows 
hydrogen at early times but then becomes helium-dominated later on (e.g., SN 1993 J; Matheson et 
al. 2000). It is possible, then, to roughly rank core-collapse events in order of increasing hydrogen 
envelope mass at the time of explosion: Ic, lb, lib, Iln, II-L and, finally, II-P. The exact sequence, 
especially the relative positions of lib, Iln, and II-L, is still quite uncertain. 

The geometry of SNe II-P are of particular interest due to their cosmological utility as extra- 
galactic distance indicators through the Expanding Photosphere Method (EPM; Kirshner & Kwan 
1974; Schmidt et al. 1994a [S94]), a variant of the Baade method used to measure distances to 
variable stars (Baade 1926). This technique is especially attractive since it is a primary distance in- 
dicator, unlike the "standard candle" method used to determine distances to SNe la. It is therefore 
not subject to all the uncertainties that plague the various rungs of the extragalactic distance ladder. 
Thus far, EPM has been applied to 18 SNe II spanning distances from 0.049 Mpc (SN 1987A) to 
180 Mpc (SN 1992am; Schmidt et al. 1994b). Using SNe II alone, a Cepheid-independent value of 
Hq = 73± 7 km s^^ Mpc^^ has been derived (S94). However, unlike the empirically-based method 
used to measure distances to SNe la, distances derived to SNe II-P rely on the assumption of a 
spherically symmetric flux distribution during the early stages of development (i.e., the plateau). 
In this way, spectropolarimetry of SNe II-P provides a critical test of the cosmological utility of 
these core-collapse events. 

Due largely to the difficulty of obtaining the required signal-to-noise (S/N) ratio, detailed 
spectropolarimetric studies exist for only three events, SN 1987A (Jeffery 1991a and references 
therein), SN 1993J (Trammell, Mines, & Wheeler 1993; Hoflich et al. 1996; Tran et al. 1997), and 
SN 1998S (Leonard et al. 2000a; Wang et al. 2000), none of which were of the classic II-P variety. 
SN 1987A resulted from the explosion of a blue supergiant, and had a peak intrinsic polarization 
of p ~ 0.7% during the first 100 days. SN 1993J, a Type lib event, was polarized at p ~ 1.5% 
at early times. SN 1998S, a peculiar Type Iln event, was polarized up to 3% at early times. The 
polarization seen in these three objects suggests that core-collapse events, at least those of the 
unusual varieties represented by them, may indeed be aspherical at early times. It is also curious 
that SN 1993 J and SN 1998S, which had both lost much of their hydrogen envelopes prior to 



exploding, were significantly more polarized than SN 1987A, which likely had retained most of its 
envelope. This trend is also seen in the sample of broadband polarimetric studies (or low S/N ratio 
spectropolarimetry) that have been carried out on SNe (Wang et al. 1996; see Wheeler [2000] for 
a comprehensive list of polarimetric observations of SNe), and tentatively suggests that the degree 
of polarization increases with decreasing envelope mass (Wheeler 2000). In fact, a polarization of 
nearly 7% has recently been inferred for a SN Ic (Wang et al. 2000a). 

It seems that the deeper we peer into the SN atmosphere, the more evidence there is for 
asphericity. This leads naturally to the conclusion that it is the collapse or explosion mechanism 
itself that is strongly asymmetric, and that evidence for this asymmetry is damped by the addition 
of outer envelope material (e.g.. Chevalier & Soker 1989). Since a SN II-P gradually reveals its 
inner layers as the recombination front recedes through the hydrogen envelope during the plateau 
phase, a single object observed at multiple epochs during (and after) the plateau era provides an 
important test of this hypothesis: explosion asymmetry should result in increasing evidence for 
asphericity with time. With the larger aperture telescopes capable of generating high S/N ratio 
spectropolarimetry in reasonable observing times now available (e.g.. Keck, VLT), the time is ripe 
for a detailed investigation into early-time SN II-P geometry. 

Supernova (SN) 1999em was discovered on 1999 October 29 UT^ by Li (1999) at an unfiltered 
magnitude of m ~ 13.5 mag in the nearly face-on [i ~ 31°, from LEDA^) SBc galaxy NGC 1637 
(Fig. 2). It was quickly identified as a Type II event, with prominent P-Cygni features of hydrogen 
Balmer and He I A5876 and a blue continuum (Jha et al. 1999). Immediately after discovery a spec- 
tral, spectropolarimetric, and photometric campaign was initiated at Lick Observatory. Nearly daily 
optical photometry and spectroscopy of the young SN were obtained with the Katzman Automatic 
Imaging Telescope (KAIT) and the Lick 1-m Nickel reflector, respectively, providing unprecedented 
temporal coverage of this Type II event during the photospheric phase. Spectropolarimetry was 
obtained using the Lick 3-m and Keck-I 10-m telescopes. 

We focus here on the spectropolarimetric observations of SN 1999em, and defer a more detailed 
discussion of the photometric and spectral data to a separate paper (Leonard et al. 2001), in which 
the EPM distance to SN 1999em (d f» 6 Mpc) is derived. For our present purposes, we note that 
both its light-curve shape and spectral properties firmly establish SN 1999em as a Type II-P event, 
with a plateau lasting roughly 90 days from discovery. In addition, the EPM analysis determines 
the time of explosion to be to ~ 5 days before discovery, consistent with the absence of the SN in 
a CCD image of NGC 1637 taken 9 days before discovery (limiting magnitude about 19; Li 1999). 
We note that at an estimated distance of ~ 6 Mpc, SN 1999em reached an extinction-corrected 
{E[B -V] = 0.05 mag [Baron et al. 2000]) peak 5-band magnitude of Mb = -15.3 ± 0.2 mag, 
which is somewhat fainter than average for SNe II-P (Patat et al. 1994). Its color evolution may 



^UT dates are used throughout this paper. 
^http://www-obs. univ-lyonl.fr/leda/home_leda. html. 



also have been somewhat unusual (Leonard et al. 2001). SNe II-P are known to be a heterogeneous 
class of objects, and these differences set a cautionary note on generalizing about the Type II-P 
class as a whole based on our findings for this one object. 

We obtained spectropolarimetry of SN 1999eni on five occasions, with observations on days 7, 
40, 49, 159, and 163 after discovery. We outline the observations and basic reduction techniques in 
§ 2 and present the data in § 3. Section 4 discusses the spectropolarimetry, the total flux spectra, 
and the implications of asphericity on the use of SNe II-P as extragalactic distance indicators 
through EPM. Conclusions are summarized in § 5. Details of the data reduction process and a 
discussion of the interpretation of low-polarization spectropolarimetry similar to what is found for 
SN 1999em are given in Appendix A. Some preliminary results of this study were announced by 
Leonard, Filippenko, & Chornock (1999) and summarized by Leonard, Filippenko, & Matheson 
(2000b). 



2. Observations and Reductions 

We obtained spectropolarimetry of SN 1999em in 1999 on November 5, December 8, and 
December 17 (7, 40, and 49 days after discovery, respectively) with the Kast double spectrograph 
(Miller & Stone 1993) with polarimeter at the Cassegrain focus of the Shane 3-m telescope at Lick 
Observatory. Similarly, on 2000 April 5 and 9 (159 and 163 days after discovery, respectively), 
we used the Low-Resolution Imaging Spectrometer (Oke et al. 1995) in polarimetry mode (Cohen 
1996^) at the Cassegrain focus of the Keck-I 10-m telescope. The slit position angle (P.A.) was 
generally set near the parallactic angle (Filippenko 1982) for the midpoint of each set of observations 
in order to minimize differential light loss. A journal of observations is given in Table 1. 

One-dimensional sky-subtracted spectra were extracted optimally (Home 1986), with a width 
of from 3 to 5 times the full width at half maximum (FWHM) of the object's spatial profile (i.e., 
the "seeing," although guiding errors and dome turbulence can also contribute). This extraction 
width was found to represent a good compromise between the increased signal-to-noise (S/N) ratio 
obtained with wide extractions and the reduction of galaxy light contamination achieved using a 
narrow aperture.^ We subtracted a linear interpolation of the median values of background windows 
on either side of the object from the object's spectrum; since the polarimetry optics restrict the 
useful (unvignetted) slit length to only 40" at Lick and 20" at Keck, background regions were 
necessarily quite narrow and close to the object. Each spectrum was then wavelength and flux 



^Instrument manual available at http://www2.keck.hawaii.edu:3636/. 

^It is also advantageous to use a large extraction window for spectropolarimetry since spurious polarization features 
can be introduced by interpolating the counts in fractional pixels at the edges of the extraction aperture. Experiments 
with different extraction widths reveal that the amplitude of such features is substantial (up to 0.5% in a Stokes 
parameter) when a width of 1 — 2 times the seeing is used, but becomes negligible (< 0.01%) when the extraction 
window exceeds 3 times the seeing. 
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Fig. 2. — S-band image of NGC 1637 taken on 1999 October 31 with the Katzman Automatic 
Imaging Telescope (Treffers et al. 1997; A. V. Fihppenko et al., in preparation). SN 1999em 
(arrow) is measured to be 15. "1 west and 17. "2 south of the galaxy nucleus (cf., Jha et al. 1999). 



TABLE 1 

Journal of Spectropolarimetric Observations of SN 1999em 



Dav" 


UT Date 


Tel.'' 


P.A.<^ 


Par. P.A."" 


Air" 


Seeing'' 


Slit 


Grating* 


Exp. 


Observer(s)'' 




(yyyy-mm-dd) 




(deg) 


(deg) 




(arcsec) 


(arcsec) 




(s) 




7 


1999-11-07 


Lick 


170 


164-184 


1.3-1.4 


2.0-2.5 


2 


600/4310, 600/5000 


4 X 1200 


AF, DL, RC 








20 


14-33 


1.3-1.6 


2.0-2.5 


2 


600/4310, 600/5000 


4 X 1500 


AF, DL, RC 


40 


1999-12-08 


Lick 


16 


22-36 


1.3-1.7 


3 


3 


600/5000 


4 X 1200 


DL, DA 








43 


42-47 


1.7-2.5 


4 


3 


600/5000 


4 X 1200 


DL, DA 


49 


1999-12-17 


Lick 


161 


156-183 


1.3-1.4 


2-2.5 


2 


600/5000 


4 X 1800 


DL, RC, AC 








24 


12-36 


1.3-1.6 


2.5 


2 


600/5000 


4 X 1800 


DL, RC, AC 


159 


2000-04-05 


K-I 


66 


64-67 


2.0-2.6 


1.6-2.1 


1.5 


600/5000 


3 X 900 


AF, AB 


163 


2000-04-09 


K-I 


66 


64-66 


1.9-2.3 


1.3-1.7 


1.5 


300/5000 


4 X 300 


BE, MB, RW 



"Days since discovery, 1999 October 29 UT (HJD 2,451,480.9), rounded to the nearest day. 
Lick = Lick 3-m/Kast Double Spectrograph with Polarimeter; K-I = Keck-I 10-m/Low Resolution Imaging Spectrometer with Polarimeter. 

"Position angle of the spectrograph slit. 

''Parallactic angle range over the course of the exposures. 

"Airmass range for each set of observations. 
The range in the FWHM of the spatial profile for each set of observations. 

^Grating denoted by (number of grooves/mm)/(blaze wavelength. A); a dichroic with crossover at 5500 A was used in the day 7 observation 
to separate the blue and red beams. 

''AB = Aaron Barth; AC = Alison Coil; AF = Alex Filippenko; BE = Bob Becker, DA = David Ardila, DL = Douglas Leonard; MB = 
Michael Brotherton; RC = Ryan Chornock; RW = Richard White. 
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calibrated, corrected for continuum atmospheric extinction and telluric absorption bands (Wade 
(fe Home 1988), and rebinned to a 2 A bin~^ linear scale. Details of the extractions are given in 
Table 2. 

We performed polarimetric analysis according to the methods outlined by Miller, Robinson, 
& Goodrich (1988) and Cohen et al. (1997). The P.A. offset between the half-wave plate and the 
sky coordinate system was determined for each observing run using polarized standard stars from 
the hsts of Hiltner (1956), Mathewson k Ford (1970), Hsu & Breger (1982), and Schmidt, Elston, 
&: Lupie (1992b). The absolute P.A. adopted for the standards is given in Table 3, along with the 
polarization angle measured for other polarized standards observed, but not used to determine the 
offset. When multiple polarization standards were observed on the same night, they always agreed 
to within ±0.5°. 

There were a few peculiarities encountered in the reductions. Briefly, on day 7, we used a 
dichroic filter that split the spectrum at about 5500 A and sent the light to separate gratings that 
dispersed the spectra in the blue and red arms of the spectrograph. This yielded spectral coverage 
from 3400 A to 7700 A, but the region around the dichroic cutoff (5100 — 6000 A) suffered from 
irregularities at the 1% level that did not calibrate out of the polarimetry; thus, we discarded the 
range 5100 — 6000 A in the polarimetric analysis. In addition, instrumental polarization ranging 
up to ~ 0.4% was detected at blue wavelengths (A < 5000 A) in this setup and had to be removed 
using the observation of the null standard on this night. On days 40 and 49 we observed without 
the dichroic and obtained continuous coverage over 4300 — 7000 A using only the red arm of the 
spectrograph. Observations of the null standards HR 7783 and HR 8086 on day 49 showed a flat, 
well-behaved polarization response. A small polarization (< 0.2%) was measured to both nulls. 
However, observations of these "nulls" on nights when additional null standards were observed 
suggest that much of this polarization is likely intrinsic to the stars themselves (i.e., the interstellar 
medium along the line-of-sight [1-o-s]) and not of instrumental origin. Therefore, we did not apply 
any instrumental correction to the data from these nights. Observations of null standards at Keck 
were always found to be null to within 0.1%, suggesting negligible instrumental polarization in the 
day 159 and 163 data. On day 159, a CCD readout failure occurred during the observation of 
the final waveplate position (67.5°), leaving us with only 3 of the desired 4 observations in the set. 
While it was still possible to derive the Stokes parameters, the data had potential systematic errors. 
Fortunately, SN 1999em was observed again just four days later at Keck, largely confirming the 
results of the earlier observation. A complete discussion of these issues, along with an investigation 
into sources of systematic uncertainty in measurements of SN polarization, is given in Appendix A. 

The displayed polarization p is the "rotated Stokes parameter" (RSP), calculated by rotat- 
ing the normalized Stokes parameters q and n by a fit to the continuum P.A. curve so that 
all of the continuum polarization falls in a single Stokes parameter (here, rotated q). This is 
preferable to calculating either the "traditional" p {ptrad = ^Jif' + "U^), or a "debiased" p {pdeb = 
i\/\ q'^ + u"^ — {Cg + o'S) I; where P(ie6 assumes the sign of [q"^ + u'^ — a'^ — a^]; see Stockman Sz Angel 



1978), since the former is biased high and is positive-definite while the latter gives poor results for 
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Table 2. REDUCTION INFORMATION FOR SN 1999EM SPECTROPOLARIMETRY 



Day'^ 


Range*^ 


Res.<= 


Flux Std. 


Ext 


. Width^ 


Background*' 




(A) 


(A) 




(arcsec) 


(arcsec) 


7 


3400-5560, 4940-7730 


5,6 


Feige 34, HD84937 




8 


-18.4:-6.4,6.4:18.4 




3400-5560, 4940-7730 


5,6 


Feige 34, HD84937 




8 


-18.4:-6.4,6.4:18.4 


40 


4230-7010 


7 


HD84937 




12.8 


-16:-8,8:16 




4230-7010 


7 


HD84937 




12.8 


-15.2:-8.8,8.8:15.2 


49 


4220-6882 


6 


HD19445 




9.6 


-20:-6.4,6.4:12 




4220-6882 


6 


HD19445 




9.6 


-20:-6.4,6.4:12 


159 


4338-6872 


7 


HD84937 




6.3 


-7.3:-4.2,4.2:5.8 


163 


3720-8685 


13 


Feige 34 




6.3 


-7.3:-4.2,4.2:5.8 



^Days since discovery, 1999 October 29 UT (HJD 2,451,480.9), rounded to the nearest day. 

'^Observed wavelength range of spectrum. In some cases, the ends are very noisy, and are 
not shown in the figures. On day 7, the spectropolarimetric data suffered irregularities near 
the crossover point (~ 5500 A) of the dichroic. The range 5100 — 6000 A is not included in 
the polarimetry analysis for this night, but is still presented in total flux since the error in flux 
calibration was only < 1% in the region. 

''Approximate spectral resolution determined by the FWHM found for night sky lines. 

'^Width of object extraction window, centered along a spline fit to the centroid of its spatial 
profile on the CCD chip. 

''The boundaries of the background windows on either side of the object, taken below and above 
the object on the CCD frame. Polarization of the diffuse host galaxy light possibly affected the 
overall measured polarization level by ~ 0.1%; see Appendix A for details. 
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Table 3. POLARIZED STANDARD STARS OBSERVED DURING NIGHTS WITH 

SN 1999EM OBSERVATIONS 





Date 






P 


P.A. 


Star 


(yyyy-mm-dd) 


Band 


Tele."^ 


(percent) 


(deg) 


HD 14433 


1999-11-05 


Sb 


Lick 


3.7 


112.9 


HD 19820 


1999-11-05 


Sb 


Lick 


4.6 


115.7 (*) 




1999-12-17 


V 


Lick 


4.8 


114.9 (*) 


HD 25443 


1999-11-05 


Sb 


Lick 


5.3 


133.9 




1999-12-08 


V 


Lick 


5.3 


134.2 (*) 


HD 43384 


2000-04-05 


V 


KI 


2.9 


170.5 (*) 


HD 204827 


1999-12-17 


V 


Lick 


5.4 


58.6 


HD 245310 


2000-04-09 


V 


KI 


4.3 


146.0 (*) 



Note. — Polarization standards used to determine the P.A. offset 
between the half-wave plate and the sky coordinate system are denoted 
with an asterisk; additional standards observed as a consistency check 
are listed here with their measured values. 

^Lick = Lick 3-m/Kast Double Spectrograph with Polarimeter; K I = 
Keck-I 10-m/Low Resolution Imaging Spectrometer with Polarimeter. 

"^Corrected for instrumental polarization; see Appendix A. 
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data with low p/op (e.g., Miller et al. 1988). RSP provides a good measure of p when the P. A. 
is a slowly varying function of wavelength. However, when there are sharp P.A. rotations across 
line features, RSP always underestimates the true polarization, and it is necessary to examine both 
RSP and URSP (rotated u) to determine the line's polarization. To measure absorption-line and 
emission-line polarizations, we first binned the flux data into 10 A (rest) bins (about 1.5 resolution 
elements wide for the standard Lick and Keck setups), and then derived q and u from the binned 
data. Although binning always reduces peak polarization, it is more reliable because of the reduced 
noise, which is especially problematic here since polarization peaks are often associated with line 
troughs. However, since features may be quite sharp, it often happens that the polarization change 
is detectable in only a few (rebinned) pixels. In such cases, we carefully examine the polarization 
properties of the original spectrum (always binned 2 A bin^^) to see if the polarization feature 
is seen consistently across many pixels or is due to a single deviant pixel, before concluding that 
a true polarization change exists. We determined observed B-hanA and F-band polarization by 
calculating the "debiased" , flux- weighted averages of q and u over the intervals 3950 — 4900 A and 
5050 — 5950 A, respectively, converting to p and 9 only at the final step for quoting results. A more 
extensive description of the adopted analysis techniques is given in Appendix A. 



3. Results 

Figures 3, 4, 5, and 6 show the observed polarization data for SN 1999em on days 7, 40, 49, and 
161 (a combination of the observations on day 159 and 163) after discovery. The data from days 7, 
40, and 49 are the weighted averages of two successive observations taken on those nights, and the 
"day 161" data are the weighted average of data obtained on days 159 and 163; see Appendix A for 
details. A montage of all the epochs studied is shown in Figure 7. Measured values for the observed 
F-band polarization are given in Table 4, along with line trough polarization for Fe II A5169 and 
Na I D AA5890, 5896, and the approximate polarization change measured across the Ha profile. 

The day 7 flux spectrum shows features that are typical for a very young SN II-P: a blue 
continuum with P-Cygni lines of hydrogen Balmer and He I A5876. The blue edges of absorption 
troughs indicate velocities in excess of 10, 000 km s^ for the fastest moving gas. The spectra from 
days 40 and 49 are typical for a SN II-P during the plateau phase, characterized by a thermal con- 
tinuum with P-Cygni hydrogen Balmer, Fe II AA4924, 5018, 5169, Sc II AA5526, 5657, 6245, and Na I 
D AA5890, 5896. The plateau phase for SN 1999em lasted until about 90 days after discovery, sug- 
gesting that a massive hydrogen envelope remained at the time of explosion (Leonard et al. 2001). 
The day 161 spectrum is emission-line dominated, with prominent Ha and [Ca II] AA7291,7324, 
implying that the transition to the nebular phase was underway; P-Cygni absorptions still indicate 
significant continuum light, however, and the fully nebular phase is not expected until after about 
day 300 (e.g., Jeffery 1991a). 



TABLE 4 
SN 1999EM: Polarization Data 



Broad-band'' 



Fe II (A5169 A) 



Na I D (AA5890, 5896 A) 



Hq 



Day* 



(percent) 



My 

(percent) 



Pv 
(percent) 



fy 

(deg) 



(percent) (percent) 



^'^^ max 

(percent) 



Vjixhr rnax 

(percent) 



(percent) 



7 


0.15 ±0.01 


-0.10 ±0.01 


0.18 ±0.01 


162.8 ±2.1 










40 


0.20 ±0.02 


-0.13 ±0.02 


0.24 ±0.02 


163.3 ±2.6 


0.4 ±0.3 


-0.1 ±0.3 


0.3 ±0.3 


0.1 ±0.3 


49 


0.27 ±0.01 


-0.25 ±0.01 


0.36 ±0.01 


158.3 ±1.0 


0.9 ±0.2 


-0.2 ±0.2 


0.5 ±0.2 


0.0 ±0.2 


161 


0.38 ±0.03 


-0.37 ±0.03 


0.53 ±0.03 


158.1 ±1.6 


1.6 ±0.7 


-3.5 ±0.6 


3.4 ±1.4 


-1.2 ±1.0 



-0.07 
0.28 
0.32 
0.56 



Note. — Quoted errors are \a statistical. Systematic uncertainties may be as high as 0.2% in a Stokes parameter; see Appendix A. 

''Days since discovery, 1999 October 29 UT (HJD 2,451,480.9), rounded to the nearest day. 

•^F-band polarization values measured by calculating the debiased (see Appendix A), flux- weighted averages of q and u over the interval 
5050-5950 A. Due to its different spectral range, the broad-band values for day 7 were taken from 6050 — 6950 A. 

•^The maximum RSP in any 10 A bin within 1000 km s~^ of the flux minimum. 

''The value of URSP for the 10 A bin in which RSP peaks. 

''The change in polarization across the Ha profile, defined as Ap = pi — P2, where pi and p2 represent the simple averages of the RSP 
values of 10 A bins from 30 to 80 A blueward and redward of the Hq flux minimum, respectively. 
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Fig. 3. — Polarization data for SN 1999em, obtained 1999 November 5, only 7 days after discovery. 
In this, and all subsequent figures, the NED recession velocity of 717 km s^^ has been removed. 
Instrumental polarization at blue wavelengths has been removed from the data on this night only 
(see Appendix A), (a) Total flux, in units of 10~^^ ergs s~^ cm~^ A~^. (6) Observed degree of 
polarization, (c) Polarization angle in the plane of the sky. (d, e) The normalized q and u Stokes 
parameters. (/) Average of the (nearly identical) \a statistical uncertainties in the Stokes q and u 
parameters for the displayed binning of 2 A bin~^. The polarization shown in this and all plots 
is actually the "rotated Stokes parameter" (RSP; see Appendix B). To facilitate comparison, the 
ordinate ranges of panels (6), (c), (d), (e), and (/) are the same in Figures 3, 4, and 5. 
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Fig. 4. — As in Figure 3, but for data obtained on 1999 December 8. 
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Fig. 5. — As in Figure 3, but for data obtained on 1999 December 17. 
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Fig. 6. — As in Figure 3, but for the weighted average of polarization data obtained on 2000 April 
5 and 9, 159 and 163 days after discovery, respectively. 
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Fig. 7. — Montage of the observed optical polarization of SN 1999em, with arbitrarily scaled total 
flux spectra {thick lines) overplotted for comparison of features. The polarization spectra are binned 
5 A bin~"^ to improve the S/N ratio. 
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4. Analysis and Discussion 
4.1. Spectropolarimetry 

4.1.1. Introduction 

We interpret the spectropolarimetry of SN 1999em in terms of tlie pliysical picture sliown 
(for a splierical atmosphere) in Figure 1, and described in § 1. Early theoretical work on SN 
spectropolarimetry was stimulated by observations of SN 1987A. Investigators initially proposed 
that the observed polarization (up to pv ~ 0.7% during the first 100 days after explosion) was 
due to the polarizing effect of line scattering (Jeffery 1989). It was soon realized, however, that 
scattering by free electrons extending out into the line-forming regions of the envelope dominates 
the overall polarization level of photospheric SNe (II91), and that weak collisions in the line- forming 
region of the SN atmosphere tend to destroy the intrinsic line polarizing effect, except perhaps in 
the outer regions of the atmosphere at late times (Jeffery 1991b). The loss of geometric information 
caused by the collisions generally produces unpolarized emission components for the P-Cygni line 
profiles, an assumption often used to constrain the interstellar polarization (ISP; see, e.g., Leonard 
et al. 2000a). Detailed models specific to SN 1987A suggest that the net polarization seen was 
likely produced by a modest asphericity of about 10% — 20% (Jeffery 1991b; H91; see also Wang &: 
Wheeler 1996). 

Translating observed SN continuum polarization into an estimate of the asphericity of the 
electron-scattering atmosphere is quite model-dependent. The original problem of the polarization 
of light emerging from a semi-infinite, plane-parallel, pure electron scattering atmosphere with 
constant net flux was solved by Chandrasekhar (1946). The first attempt to apply the formalism 
to SNe was by Shapiro & Sutherland (1982), who solved it for a plane-parallel, optically-thick, 
electron-scattering atmosphere wrapped around an ellipsoid. The plane-parallel model, however, 
requires extreme asymmetry to produce even modest polarization: to get p = 1% demands an 
asphericity of at least 60% (the percent asphericity is defined here as (1 — [a/b]) x 100, where a 
and b are the semiminor and semimajor axis lengths, respectively) for the oblate case; prolate 
spheroids are unable to produce polarization of more than 0.77% for any axis ratio or viewing 
angle. Realistic models of the extended atmospheres of SNe, however, have recently been shown to 
be more polarizing than the plane-parallel approximation predicts. 

There are many factors that can influence the degree of SN continuum polarization. The most 
recent, general study of the polarization expected from an aspherical SN atmosphere is H91, in which 
the continuum polarization resulting from an axially symmetric (either oblate or prolate), electron- 
scattering dominated photosphere is calculated as a function of the photosphere's geometry and size 
relative to the envelope, the scattering atmosphere's optical depth and density profile, and the axial 
ratio of the envelope. Particularly useful is the lower bound placed on the asphericity implied by 
polarization measurements: by fixing all other parameters to maximize the resulting polarization 
(i.e., scattering optical depth = 1, continuum region physically small compared with envelope. 
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object viewed equator-on), the maximum polarization as a function of axis ratio is produced (Figure 
4 of H91). Determination of intrinsic SN continuum polarization can thus lead to a lower limit on 
the implied asphericity of the scattering atmosphere from this model. 

In addition to the information gleaned from the continuum polarization level, sharp polariza- 
tion changes across strong line features can also contain a wealth of information about the geometry 
of the scattering environment. The influence of asphericity on the polarization properties of P-Cygni 
lines was first addressed in detail by McCall (1984), who predicted that polarization should increase 
in the absorption troughs. His basic idea is that P-Cygni absorption by gas in region 1 (Fig. 1) 
selectively blocks photons coming from the central, more forward-scattered (hence less polarized) 
regions, thereby enhancing the relative contribution of the more highly polarized photons from the 
limb (i.e., from regions 2 and 3 in Fig. 1; see also Fig. 1 of Leonard et al. 2000b). Jeffery (1991b) 
extended the argument by predicting drops in polarization at the location of emission peaks due 
to the dilution of the polarized continuum light by unpolarized emission-line photons. 

The basic prediction, then, is that a "reverse" P-Cygni polarization profile should characterize 
strong P-Cygni line features formed above or within an aspherical electron-scattering atmosphere. 
(We note that polarization increases in P-Cygni absorption troughs are often seen in spectropolari- 
metric studies of broad absorption line QSOs, and are generally given a similar physical interpre- 
tation [e.g., Brotherton et al. 1997; Ogle et al. 1999]). Taking the limiting case of blockage of only 
unpolarized light (or light with electric vectors that exactly cancel) in a line trough and addition 
of completely unpolarized light in a line peak, the polarization change expected across a P-Cygni 
feature from this admittedly simple model is 

^ ^cont /-. \ 

P trough — Prnax — Pcont ^ T i {^) 

J^ trough 

^ ^COnt /r\\ 

Ppeak ^ Pmin — Pcont ^ "7 ) \^) 

^peak 

where Pcont is the continuum polarization and Icont, Ipeak, and hrough represent the flux intensity 
in the continuum, line peak, and line trough, respectively. The continuum polarization and the 
continuum flux intensity are both interpolated from the continuum regions surrounding the line. 

In general, the polarization observed across a P-Cygni profile will differ from that predicted by 
equations (1) and (2). The most important factor neglected by the simple geometric argument is 
the influence that resonance scattering of continuum photons has on the polarization properties of 
the line profile. In a SN atmosphere, continuum photons that are resonantly scattered by a line are 
unpolarized (Jeffery 1991b; Hoflich et al. 1996), and will diminish the polarization level throughout 
the entire P-Cygni profile. This depolarizing effect will be especially pronounced in lines that are 
optically thick throughout the envelope, such as Ha. In addition, the strong emission component of 
Ha will tend to fill in the absorption trough, increasing Itrough- Countering this is the polarization 
that line photons themselves may acquire upon being scattered by free electrons, which likely exist in 
the line- formation region itself (Hoflich et al. 1996 ). A flnal reason that the polarization measured 
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will differ from that predicted is purely observational: since the P-Cygni absorption features in 
SNe II-P are quite sharp (the FWHM of the absorptions in SN 1999em due to prominent metal 
lines on day 49, for example, are all in the range 35 — 45 A), the typical grating-CCD resolution 
of 5 — 10 A coupled with the coarse binning needed to improve the statistics inevitably reduces 
measured peak polarization values. Indeed, in the limit of infinite bin size, Ptrough -^ Pcont- 

There are thus several mechanisms in addition to the simple geometric model that may con- 
tribute to the polarization measured across a P-Cygni line feature and it is not obvious which 
should dominate. In this regard, it is instructive to study the characteristics of the line polar- 
ization presented by the detailed spectropolarimetric studies of SN 1987A (Jeffery 1991a, 1991b), 
SN 1993 J (TrammeU et al. 1993; Hoflich et al. 1996 ; Tran et al. 1997), and SN 1998S (Leonard 
et al. 2000a). 

For SN 1987A, many epochs of the ISP-subtracted data show large polarization increases at 
the locations of P-Cygni absorptions due to H7, B.f3, Fe II A5169, Na I D AA5890, 5896, Ca II 
AA8498, 8542, 8662, and possibly Fe II A5018 and Ha as well (see Jeffery 1991a and references 
therein). Polarization decreases occur across the emission peaks of Ha and Ca II AA8498, 8542, 8662. 
Since error bars are not given and the polarization data were evidently not corrected for statistical 
biases, it is somewhat difficult to quantitatively assess the polarization changes (especially the 
trough increases) seen across the lines. Taken at face value, though, none of the trough increases 
appears to exceed the limit given by equation (1). Polarization decreases across Ha, however, do 
occasionally go below the limit set by equation (2). 

The three ISP-subtracted spectropolarimetric epochs of SN 1993 J (Type lib) presented by Tran 
et al. (1997) also show distinct polarization increases in the absorption troughs of H/3, Fe II A5169, 
and He I A5876 (probably blended with Na I D AA5890, 5896), and polarization decreases across 
the He I A5876 and blended Ha/He I A6678 emission lines (the reverse P-Cygni profile across the 
He I A5876 is particularly pronounced). Detailed modeling of a single early-time spectropolarimetric 
epoch of SN 1993J presented by Hoflich et al. 1996 successfully reproduces the depolarization across 
the entire Ha profile, but does somewhat less well at fitting the He I/Na I D feature: whereas the 
model predicts a polarization decrease at the line's absorption trough, the data indicate a small 
polarization rise. It is possible that the increased importance selective blocking of forward-scattered 
light has in lines of low optical depth contributes to the different polarization behavior seen here 
between Ha and He I/Na I D. 

Perhaps the clearest demonstration of the influence that unpolarized emission lines have on 
SN spectropolarimetry, however, is the single spectropolarimetric epoch of SN 1998S discussed by 
Leonard et al. (2000a). At the very early epoch studied (just 5 days after discovery), this peculiar 
Type Iln event displayed a smooth, blue continuum with superposed broad, symmetric emission 
lines lacking any absorption components. The spectropolarimetry was characterized by a fiat, 
highly polarized (p ~ 3% after ISP removal) continuum with pronounced polarization decreases 
across all the emission features. 
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Aspects of the predicted reverse P-Cygni polarization profile are thus found in the spectropo- 
larimetry of all previously studied SNe II. However, SN 1987A, SN 1993J, and SN 1998S were 
all quite unusual core-collapse events. SN 1993J and SN 1998S both showed strong spectroscopic 
evidence for interaction with circumstellar material. In addition to being a blue supergiant, the 
progenitor of SN 1987A may also have been spun up by a merger with a companion star prior 
to explosion (Collins et al. 1999). Such circumstances add complexity to the interpretation of 
spectropolarimetry. Since SN 1999em appears to be a relatively normal SN II-P, it presents a great 
opportunity to study the impact of asphericity on line polarization without any added complication 
by "external" influences. 



4.1.2. Interstellar Polarization 

A major source of uncertainty in SN spectropolarimetry is the degree to which interstellar dust 
along the 1-o-s polarizes the SN light. An assumption often invoked to determine the ISP to SNe 
is that broad, resonance-scattered and recombination line photons (e.g., Ha) at early phases are 
intrinsically unpolarized; see § 4.1.1 for the theoretical underpinnings behind this assumption, and 
Appendix B for details of the application of this technique (see also Jeffery 1991b; Trammell et al. 
1993; Hoflich et al. 1996 ; Wang et al. 1996; Tran et al. 1997). With this assumption, the observed 
polarization of the emission-line photons is then entirely due to the ISP. Our earliest epoch, day 
7, has very low continuum polarization and no obvious line features (see Fig. 12 for detail of the 
Ha region). For Ha emission from 6404 to 6622 A (± 5000 km s^ from the deredshifted line 
peak at 6512 A) and background regions 6000 — 6250 A and 6900 — 7500 A, we find p{Ha) = 
0.08 lb 0.06% at ^ = 171° ± 15°. This value is consistent with the continuum measured in the range 
6900 - 7500 A (p = 0.12 ± 0.02% at 6^ = 153 ± 5°), and therefore allows for the possibility that all 
of the polarization seen is interstellar in origin. On the other hand, the polarization found is also 
marginally consistent with zero polarization, which would allow all of the observed polarization to 
be intrinsic to the SN. Our data therefore do not distinguish well between an ISP of 0.12% and an 
ISP of zero with this test. In any event, though, if Ha emission-line photons are really intrinsically 
unpolarized at this early epoch, it does suggest very low ISP from dust in either the host galaxy 
or Milky Way (MW), with a formal 3a upper hmit of ISP < O.z 



There are additional indications of low ISP contamination along the 1-o-s to SN 1999em. First, 
since the same dust that polarizes also reddens SN light, with the maximum polarization always 
found to be less than 9E{B - V) for the MW (Serkowski, Mathewson, & Ford 1975), a small ISP 
for SN 1999em is consistent with its low Galactic reddening of E{B — V)mw = 0.04 mag (Schlegel, 
Finkbeiner, & Davis 1998), as well as with the total reddening of E{B — V)tot ~ 0.05 mag estimated 
by Baron et al. (2000; see also Leonard et al. 2001) from theoretical fits to the early-time continuum 
shape. This reddening restricts the ISP to values less than 0.45%; we note, however, that the upper 
limit of E(B—V) = 0.15 mag given by Baron et al. (2000) is consistent with an ISP of up to 1.35%. 
Second, of the 11 stars with measured polarizations within 5° of SN 1999em (Heiles 2000), six have 
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zero polarization within their quoted uncertainties, and all have pv < 0.2%, although only two 
(HD 29248 and HD 30112) he sufficiently far away {d > 300 pc ) to fully sample the ISP through 
the Galactic plane. Third, even though p increases through the 4 epochs, the polarization angle 
of SN 1999em remains nearly constant at ^ ~ 160°. If the day 7 polarization was completely due 
to ISP it would suggest that Ofsp = ^SNi an unlikely chance occurrence. Furthermore, optical 
imaging polarimetry of face-on spirals (Scarrott, Rolph, & Semple 1990; Scarrott et al. 1991) 
generally shows polarization vectors aligned parallel to the spiral arms (i.e., perpendicular to the 
line connecting a point with the galaxy's center), which would predict 9isp ~ 130° given the 
location of SN 1999em in NGC 1637, somewhat far from the observed 9 ~ 171°. Finally, the nearly 
complete depolarization across the Ha line to q ^ 0%,u « 0% on days 40, 49, and 161 (see Figs. 
13, 14, 15) is consistent with a near-zero ISP from physical considerations discussed in the next 
section. 

Therefore, while none of the preceding arguments is individually conclusive, taken together 
they suggest that polarimetry of SN 1999em is contaminated by a relatively small degree of ISP. 
Although impossible to know for sure, we suspect that the ISP is less than 0.3%, and cautiously 
proceed to analyze our data with no correction for its effects. We note that since ISP is a smoothly 
varying function of wavelength and is constant with time, uncertainty in its value affects neither 
temporal polarization changes nor modulations across specific line features. We caution, however, 
that any interpretation of polarization increases, decreases, or directional changes in line features is 
affected by the choice of ISP (see, e.g.. Fig. 10 of Leonard et al. 2000a). Other sources of systematic 
uncertainty may add an additional 0.1% uncertainty to the accuracy of the measured polarizations 
(see Appendix A). 



4.1.3. SN 1999em 

The broadband polarization of SN 1999em (generally, observed V) rises fromp = 0.18ib0.01% 
on day 7 to p = 0.53 ±0.03% on day 161 (Table 4 and Fig. 8a; quoted uncertainties are la statistical 
only, and do not reflect systematic effects) , suggesting a modest amount of intrinsic SN continuum 
polarization. The lack of any obvious wavelength dependence of the continuum polarization (e.g.. 
Fig. 7) is consistent with an electron-scattering origin. A remarkable feature is that although the 
polarization level changes, it maintains a nearly constant polarization angle in the plane of the 
sky (Fig. 8b), similar to what was seen in SN 1987A (Jeffery 1991a). This strongly suggests that 
a preferred symmetry axis persists as the photosphere recedes and exposes ever deeper regions of 
the envelope. The change in polarization by ~ 0.3% from day 7 to day 161 cannot be produced 
by ISP, and is likely intrinsic to the SN itself. We note that the overall polarization level and 
characteristics of the line features are consistent with those reported by Wang et al. (2000b) for 
spectropolarimetry of SN 1999em obtained on day 72. The continuum polarization level implies 
a minimum asphericity of ~ 7% during the plateau era and ~ 10% on day 161 from the oblate, 
electron-scattering models of H91. 
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Fig. 8. — (a) Observed broadband polarization and (6) polarization position angle in the plane of 
the sky of SN 1999em as a function of time. Error bars are la statistical for the observed V band 
(5050 — 5950 A) except for day 7, which is for observed 6050 — 6950 A. Note that although the 
polarization more than doubles during the first 161 days, the polarization angle remains fixed. 
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The temporal polarization increase may indicate increasing asphericity deeper into the SN 
ejecta. However, it is not clear that the models investigated by H91 are directly applicable for 
either very early or very late SN epochs. At early times (i.e., prior to the appearance of line 
features due to heavy metals; for SN 1999em this occurred on about day 9 (Leonard et al. 2001), 
the electron density profile above a SN II photosphere is believed to be very steep (see, e.g., Jeffery 
1991b; Eastman et al. 1996; Matzner & McKee 1999). This type of atmosphere has much less 
polarizing power than one with the shallower electron density gradients typical of SN atmospheres 
during the recombination phase studied by II91. In fact, Jeffery (1991b) suggests that the shift from 
a steep to a shallow electron density gradient may explain the polarization rise seen in SN 1987A 
early on (from pv ~ 0.2% to pv ~ 0.7% over the first 30 days) without the need for any increase 
in asphericity. Hoflich et al. 1996 cautiously propose the same mechanism to explain the early- 
time polarimetric behavior of SN 1993J as well. It is possible that a similar effect is at work in 
SN 1999em. The characteristics of the electron-scattering atmosphere during the transition to the 
nebular phase (i.e., day 161) also differ significantly from those modeled by II91. As SN 1999em 
ages the optical depth to electron scattering decreases; at some point in the nebular phase, in fact, it 
will drop well below 1, rendering spectropolarimetry ineffective as a geometric probe. It is possible 
that a reduced electron-scattering optical depth is responsible for the relatively low polarization 
measured at this late epoch. 

While interpretation of the overall polarization level is hampered somewhat by possible ISP 
contamination and other sources of systematic uncertainty (Appendix A) , the distinct polarization 
modulations across strong P-Cygni features are not. As anticipated (§ 4.1.1), sharp polarization 
peaks are associated in the later epochs with deep P-Cygni line troughs, most conspicuously Fe II 
A5169 and Na I D AA5890, 5896 (probably blended with Ba II A5853 and He I A5876 on day 161). 
We measure trough polarization in 10 A bins to help improve the statistics in these low S/N ratio 
regions, and display the results in Figures 9, 10, 11 and in Table 4. The most dramatic increase is 
seen in the blended Na I D line on day 161, in which the measured RSP is p = 3.4% it 1.4%, and 
the true polarization is likely even higher given the evidence for P.A. rotation found in its URSP 
(see Appendix B). The measured polarization increases in all troughs are in accord with the limit 
set by equation (1). 

The strong depolarization across the Ha profile also implies intrinsic SN polarization (Figs. 12 
through 15), with the change in polarization across the line increasing by ^ 0.6% during the first 
161 days (Table 4). On days 40, 49, and 161, nearly complete depolarization exists from slightly 
redward of the flux absorption minimum to the line peak near t; ~ km s^^. The fact that the 
depolarization goes well below the value expected from the simple model of dilution of polarized 
continuum light with unpolarized Ha photons (eq. [2]) and the lack of a polarization increase 
associated with the Ha trough may be due to a depressed underlying continuum polarization 
resulting from resonance scattering of continuum photons by Ha (§ 4.1.1). In effect, since Ha 
is optically thick throughout the atmosphere, it presents an absorbing screen that "covers" both 
the limb and the central regions (Fig. 1), unlike the optically thin metal lines (i.e., Fe II or Na I 
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Fig. 9. — Detail of the region around the absorption trough of Fe II A5169 on 1999 December 17. 
(a) Normalized flux, displayed at 2 A bin^^. (6) Observed degree of polarization, measured as the 
rotated q parameter, (c) Observed degree of polarization in the rotated u parameter. Error bars in 
(6) and (c) are la statistical for 10 A bin~^. In the absence of rapid P. A. change with wavelength, 
the URSP should scatter about zero; see Appendix B for details. 



28 




-10000 



-5000 5000 

Velocity (km s~^) 



10000 



Fig. 10. — As in Figure 9, but for the weighted average of data obtained on 2000 April 5 and 9. 
The sharp spike in URSP at the bottom of the hue trough may indicate that RSP underestimates 
the true trough polarization. 
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Fig. 11. — As in Figure 9, but for the weighted average of data obtained on 2000 April 5 and 9 
near the absorption trough of Na I D AA5890, 5896 (Ba II A5853 and He I A5876 probably also 
contribute). Note that URSP ~ —1.7% in many consecutive pixels near the absorption minimum, 
indicating that RSP may underestimate p by as much as 0.4%. 
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D) that predominantly block unpolarized (i.e., forward-scattered) light (§ 4.1.1). A similar Hq 
depolarization was observed and theoretically modeled in SN 1993 J (Hoflich et al. 1996 ), and was 
seen in SN 1987A as well (Jeffery 1991a; see also § 4.1.1). 

On day 161 the polarization increases steadily through the red wing of the Ha profile, reaching 
p ~ 1.5% at f = 5, 000 km s"^, and then drops back to the continuum value at w > 12, 000 km s~^. 
The polarized flux (shown overplotted in Fig. 15a) also peaks about 3200 km s~^ to the red of the 
nominal line center. One agent that may contribute is the redshift that line photons acquire upon 
being scattered by electrons moving away from the point of origin (i.e., the expanding universe 
paradigm; see, e.g., Bailey 1988; Witteborn et al. 1989). Another possibility is that at this 
late epoch the density and temperature conditions in the outer atmosphere are such that the 
depolarizing effect of weak collisions is reduced, allowing the intrinsic line polarizing effect due to 
resonance scattering to contribute (see Jeffery 1989, 1990 for a thorough analysis of polarized line 
radiative transfer in an expanding atmosphere). A final, intriguing possibility is that Ha photons 
on the "far" side of the expanding SN (region 4 in Fig. 1), are scattered into the 1-o-s by electrons 
in region 3. Although at this late epoch the physical size of region 4 is quite small and would 
contribute little flux, the photons would be highly polarized due to the near 90° scatter they suffer, 
perhaps allowing them to significantly affect the polarization of the line's red wing. 

The existence of a broad, strong, polarized emission line allows us to investigate a key question 
about the line-scattering region's geometry: is it axisymmetric? That is, does the line profile 
present a uniform polarization angle in the plane of the sky? The best way to test this is to follow 
the path traced by the data in the q-u plane, since a common axis of symmetry for the different 
scattering regions will result in a linear path. This method is preferred over simply looking to see 
if the polarization angle changes across the line (i.e., 6 in Fig. 6), since ISP can produce apparent 
P.A. rotations, but will not alter the shape of the path in q-u. Figure 16 shows the results for Ha 
on day 161, extending from to 5000 km s~^. The path traced by the red side of the profile is 
quite linear, and the polarization angle is consistent with that measured for the continuum: for line 
photons in the range v = — 5000 km s^^, we measure p = 0.48% at ^ = 161°. This measurement is 
complicated by the lack of a clear continuum, however, which required us to estimate the continuum 
level by hand in the Stokes parameter fluxes (see Appendix B). This implies that the scattering 
region responsible for the polarization exhibited by the red wing of Ha is axisymmetric and has the 
same geometry as the continuum scattering region, consistent with the idea that both the line and 
the continuum are scattered by the same electrons in the line-forming region of the SN envelope. 
Unfortunately, the low S/N ratio for v > 5000 km s~^ prohibits an investigation out to the highest 
velocities; similarly low S/N ratio in other emission lines (and Ha at other epochs) prohibits detailed 
analysis in them as well. On the blue side of the line profile, it is possible that the line prefers a 
path different from that seen on the red side, but the limited number of points preclude a definitive 
conclusion about the breaking of axisymmetry there. In all, though, the near-zero polarization of 
the Ha emission profile seen in Figure 16 near v = km s^^ bolsters the belief that the ISP may, 
in fact, be quite small. 
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Fig. 12. — Region around Ha on 1999 November 5. Error bars are la statistical for 10 A bin~^. (a) 
Total flux, normalized by a spline fit to the continuum, displayed at 2 A bin^^ for better resolution. 
(6) Observed degree of polarization, determined by the rotated Stokes parameter (rotated q). (c) 
Observed degree of polarization in the rotated u parameter, (d, e) Normalized q and u Stokes 
parameters. To facilitate comparison, the same ordinate ranges are used in Figures 12, 13, and 14. 
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Fig. 13. — As in Figure 12, but for data obtained on 1999 December 8. Note the depolarization 
evident across the emission profile. 
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Fig. 14. — As in Figure 12, but for data obtained on 1999 December 17. The depolarization seen 
across Ha on day 40 (Figure 13) is evident here as weh. 
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Fig. 15. — As in Figure 12, but for the weighted average of data obtained on 2000 April 5 and 9, and 
shown with an extended velocity range. The normalized polarized flux (RSP xfx) is overplotted 
(thin line) in (a). Note the increased polarization in the red wing of the line. 
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Fig. 16. — Polarization of the Ha line profile in the q-u plane from v = —1000 to v = 5000 km 
s^^. Symbol size increases with velocity. Error bars are la statistical for 10 A bin^^. The linear 
path traced from w = to w = 5000 km s^^ suggests an axisymmetric scattering environment for 
those line photons; the different q-u location favored by the photons blueward of v = km s~^ may 
indicate a breakdown of axisymmetry across the line profile, but the small number of points makes 
this difficult to confirm. 
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4-1 .4- The Effect of Asphericity on the Expanding Photosphere Method 

If the geometry of a SN II-P photosphere resembles the geometry of its electron-scattering 
region, then asphericity will affect the distance derived by any method that assumes a spherical 
flux source, such as EPM. Wagoner (1991 [W91]) derives the error in EPM distances resulting 
from an oblate SN photosphere with constant flux over its surface for all viewing orientations and 
degrees of asphericity. For instance, a 10% asphericity is shown to produce an EPM distance that 
overestimates the actual distance by ~ 5% for an edge-on view and underestimates it by ~ 10% for a 
face-on 1-o-s. The error naturally rises with increasing asphericity, reaching over 50% for axis ratios 
of 2:1 viewed face-on. While individual distances may thus be quite inaccurate for aspherical flux 
distributions, an encouraging result is that the average distance derived to SNe II-P with random 
viewing orientations is quite robust, with cIepm = dactual to within 2% (in the sense dEPM < dactual) 
for mean axis ratios of up to 2:1. The main effect of asphericity on the determination of Hq from 
SNe II is to increase the scatter of EPM distances about the true values, but not to substantially 
bias the overall measurement of Hq. The amount of increased scatter expected for a given mean 
SN II asphericity is shown in Figure 4 of W91. 

A useful characteristic of EPM is that each observational epoch provides an independent es- 
timate of the SN's distance once the time of explosion is derived (see, e.g., Schmidt, Kirshner, & 
Eastman 1992a). If the asphericity of SN 1999em really does increase with time while maintaining 
a constant viewing orientation (as possibly suggested by its increasing polarization at a constant 
P.A.), then the EPM-derived distance should deviate more and more from the correct value as 
time goes on. The EPM technique itself can thus be used to answer some questions concerning the 
nature of a particular SN's asphericity: EPM-derived distances that increase with time would indi- 
cate a more edge-on viewing orientation, while those that decrease with time would indicate a more 
face-on viewing orientation, with the transition from "more edge-on" to "more face-on" occurring 
at ^ 38° for the simple model of W91. (EPM will produce a non- varying distance to either a 38° 
viewing orientation or a constant degree of asphericity; subtle combinations of temporally varying 
viewing orientation and asphericity could also conspire to produce this result.) 

In this regard, it is useful to reexamine EPM distances previously derived for the 16 SNe II 
tabulated by S94, which include 10 SNe II-P and 6 SNe II-L. None of the SNe II studied shows 
any obvious temporal trends in its derived distance (e.g.. Fig. 8 of Schmidt et al. 1992a; Fig. 6 
of Schmidt et al. 1994b), with several having data points that extend to more than 60 days after 
shock breakout. In the EPM analysis of SN 1999em (Leonard et al. 2001) the derived distance 
to SN 1999em also remains quite stable with time. Barring a conspiracy of offsetting errors,^ the 
temporal constancy of EPM distances implies that any asphericity that exists must not change 



^We note that the main source of systematic uncertainty in EPM distances comes from the "dilution factor" 
used to correct the observed flux due to the departure of the photospheric flux from that expected for a blackbody 
radiating at the location of the optical photosphere. For more details, see Eastman et al. 1996, Hoflich 1991b, and 
Baron et al. 1995. 
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appreciably as a particular SN evolves during the photospheric phase. 

If SNe II are aspherical but present a fixed degree of asphericity during the photospheric phase, 
a systematic increase in the scatter in a Hubble diagram beyond that which can be explained by 
statistical uncertainties alone will result. Figure 4 of W91 shows the anticipated degree of this 
systematic scatter (given as AQrms, where Q = dA/dE and dA is the "Hubble distance" [defined 
below], dE is the EPM-derived SN distance, and AQ = Q— 1.0) as a function of mean SN asphericity. 
We may therefore gain some insight into mean SN II asphericity during the photospheric phase by 
examining the scatter of Q for the 16 SNe II reported by S94. To do this, we form a x^ statistic, 
where 

2 _ V- {Q^-l■0)^ 



X 



1 Wsys,i ^stat,i 



1 ^2 I f'^dA'A'^ I / '^A,i'ydE, 



c sys,i \ d'E , ' 



(3) 



where aq^y^ is the systematic uncertainty in Q (presumably due to asphericity), ctq^^^^ is the statis- 
tical uncertainty in Q, dE and ad^ are the "EPM distance" and statistical uncertainty, respectively, 
taken from column 6 of Table 6 by S94, dA is the "Hubble distance" {dA = Vrec/Ho, where Vrec is 
the recession velocity of the host galaxy, taken from column 3 of Table 6 by S94 and Hq is chosen 
so that the weighted mean of Q is 1.0), and a^i^ is the uncertainty in the Hubble distance due to 
peculiar velocities of the host galaxies. Peculiar galaxy velocities are generally thought to be in the 
range 100 < a^^^^ < 400 km s~^ (Lin et al. 1996; Marzke et al. 1995). We stress that the point of 
our investigation here is not to derive -f^O) but rather to set a uniform fiducial scale from which to 
quantify the scatter of EPM-derived distances about their "actual" distances. The particular value 
of Hq required to make Q = 1.0 is sensitive to both dE and dA of the SN sample, as well as the 
errors assigned to dE, dA, and Qsys- The main assumption made by setting dA = Vrec/HQ is that 
the Hubble law is linear out to 180 Mpc, the farthest EPM distance yet reported (S94; Schmidt et 
al. 1994b). 

The resulting scatter of Q about unity is shown in Figure 17 for all 16 SNe II. The la errors 
on Q and aq incorporate only the statistical uncertainty reported for dE by S94. For the full set of 
16 SNe II we derive a reduced x^ of 1.20, implying that the uncertainty resulting from ad^ alone 
will produce at least this amount of scatter 26% of the time. Including a small contribution to 
the uncertainty in Q from peculiar velocities of cTvrec — 100 km s^^ increases this likelihood to 
70%. In other words, assuming that the only sources of uncertainty in Q are the reported errors 
on dE and a small contribution from peculiar velocities leads to the conclusion that the observed 
scatter is consistent with no directional dependency of the SN fiux at the la level. If we follow 
S94 and assign a 300 km s^^ uncertainty to Vrec due to peculiar velocities, then the likelihood that 
the observed degree of scatter can be explained with no additional sources of uncertainty rises to 
99%, and the formal 3a upper limit on AQums due to directionally dependent flux is found to 
be AQrms < 0.16. Using Figure 4 of W91, this translates into a 2>a upper limit on mean SN 
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asphericity during the photospheric phase of about 30%. 

Therefore, unless the uncertainties reported by S94 do not accurately reflect (i.e., they overesti- 
mate) the true statistical error on (Ie or a conspiracy of offsetting errors exists, significant directional 
dependence of SN II luminosity is inconsistent with the existing EPM data; the subsets of SNe II-P 
and SNe II-L alone support a similar conclusion for each sample individually, although the reduced 
sample sizes produce somewhat less robust constraints. Although additional spectropolarimetric 
data are needed to create a sample large enough to directly estimate average SN II asphericity, the 
apparent success of EPM thus already significantly limits its value. The good agreement among 
EPM, Tully-Fisher, and Cepheid distances (S94; Eastman et al. 1996) also argues for low SN II 
asphericity during the photospheric phase. 



4.2. Flux Spectra 

A nice byproduct of the long exposure time required for spectropolarimetry (typically about 
4 hours to achieve (Tq^u ~ 0.3% with 2 A bin~^ using the Lick 3-m reflector while SN 1999em was 
on the plateau with ray < 14 mag) is the very high S/N ratio total flux spectrum that results. 
Although the spectral evolution of SN 1999em will be discussed in more detail elsewhere (Leonard 
et al. 2001), we comment briefly on two of the more interesting features seen in the flux spectra 
presented here. 

Figure 18 shows the absorption troughs for the He I A5876 and hydrogen Balmer lines. A 
curious "dip" exists at v ~ —12, 600 km s^^ in the higher-order Balmer lines, and possibly He I as 
well. One possibility is that these dips are produced by Chugai's (1991) proposed mechanism for 
the anomalous blueshifted "absorption" seen in SN 1987A (i.e., the "Bochum event"). Under this 
model, the features are actually produced by a lack of absorption by gas at velocities surrounding 
the "dip," caused by ineffective screening of continuum photons due to excitation stratification in 
the envelope. On the other hand, the features may be related to the "complicated" P-Cygni profiles 
reported by Baron et al. (2000) from analysis of a spectrum of SN 1999em taken within a day of 
discovery. From detailed spectral modeling, Baron et al. (2000) conclude that broad absorption 
shelves seen in the H/3 and He I A5876 profiles at u ~ 20, 000 km s~^ are due to the combination of 
a shallow density gradient and high photospheric temperature at this early epoch. Similar features 
are seen in very early spectra of SN 1990E and were identified with N II AA4623, 5679 (Schmidt et 
al. 1993), but Baron et al. (2000) favor the "complicated" P-Cygni explanation from their model 
fits. The consistent velocity of the dips in our spectrum and its presence in H7 also argue against 
the N II identification for these features. Whether the dips are in fact related to the earlier features 
seen by Baron et al. (2000) will have to await further modeling. 

Another interesting aspect is the development of a distinct "notch" centered near the rest 
wavelengths of prominent emission lines on day 161 (Fig. 19), creating a subtle "double peak" to 
the emission profiles of Hq, H/3, and [Ca II] AA7291,7324. Inspection of earlier spectra indicates 
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Fig. 17. — Ratio (Q) of the distance determined by assuming a linear Hubble law {cIa = Vrec/Ho, 
where Vrec is the host galaxy's recession velocity) to the distance determined by EPM (ds) for 
16 SNe II, vs. d^- The displayed la uncertainties in Q are derived entirely from the statistical 
uncertainty reported for cIe- To determine Q, we chose Hq so that the weighted mean of Q is equal 
to 1; for these 16 SNe II with the displayed la uncertainties this required setting Hq = 64 km s~^ 
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Fig. 18. — P-Cygni profiles observed on 1999 November 5. A dip near v ~ —12, 600 km s^ is seen 
in the higher-order Balmer hnes and possibly He I A5876 as well. Due to telluric absorption near 
6290 A (see, e.g., Matheson et al. 2000), its presence could not be verified in Ha. 



41 



1.4 



1.2 

-^^ 
a 
a 

^ 1.0 

o 
U 

+ 0.8 

^ 0.6 

u 



0.4 



0.2 



0.0 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



SN1999em: . 
Late-Time 
Line Profiles 




[Ca II] 

7291/7324 



J I I I I I I I I I I I I I I I I i_ 



-10000 10000 

Velocity (km s"^) 



Fig. 19. — Emission-line profiles at late times for SN 1999em, taken from spectra obtained on days 
159 (Ha and H/3) and 163 ([Call] AA7291,7324; Ao = 7291). The "notch" near zero velocity first 
appeared between days 95 and 124 and persists through the most recent spectrum obtained on day 
333 as well. 
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that this feature developed sometime between days 95 and 124, coinciding with the time SN 1999em 
fell off the plateau. It exists in a spectrum obtained 333 days after discovery as well. A mech- 
anism traditionally invoked to explain double-peaked profiles in SNe is interaction of the ejecta 
with a circumstellar disk of material, perhaps ejected by the progenitor prior to explosion (e.g., 
SN 1998S; see Leonard et al. 2000a). Since there are few other indications of circumstellar inter- 
action, however, we are quite skeptical of this explanation. It is also possible that slowly moving 
circumstellar material exterior to the expanding ejecta would produce an "absorption" feature near 
Ao . Alternatively, the appearance of the features at the same moment the hydrogen envelope fully 
recombined (i.e., the end of the plateau) coupled with their persistence through at least the first 
part of the nebular stage might suggest an origin within the expanding ejecta itself, perhaps due 
to an asymmetry in the line-emitting region. 



5. Conclusions 

We present spectropolarimetry of SN 1999em on days 7, 40, 49, 159, and 163 days after 
discovery, with the last two epochs combined to form a single "day 161" epoch. Our main results 
are as follows. 

1. The observed broadband polarization rises from p ~ 0.2% on day 7 to p ~ 0.5% on day 161. 

2. The continuum polarization does not exhibit any noticeable variation with wavelength. 

3. The polarization angle is ~ 160° throughout the first 161 days. 

4. From several lines of evidence, the ISP along the 1-o-s to SN 1999em is thought to be quite 
small, although values as high as 0.3% can not be ruled out. 

5. Increases in polarization of up to 0.54 it 0.2% during the plateau and 2.9 it 1.5% after the 
plateau are associated with P-Cygni absorption troughs of metal lines. 

6. Nearly complete depolarization across the Ha profile is seen on days 40, 49, and 161; the 
depolarization is neither confirmed nor rejected in the data from day 7. 

7. On day 161, the red wing of the Ha profile is polarized up to 0.7% higher than the surrounding 
continuum. The scattering medium responsible for the polarization of the red wing of the Ha 
profile traces a linear path in the q-u plane, with a polarization angle consistent with that of 
the continuum. 

8. Analysis of EPM-derived distances to 16 SNe II reveals little evidence for directional de- 
pendency of SN flux. Making reasonable assumptions about the extent of other sources of 
uncertainty in EPM distances, we derive a 3a upper limit of 30% on mean SNe II asphericity 
during the photospheric phase. 
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9. In the total flux spectra, a "dip" exists in the absorption profiles of H/3 and H7 at w ~ 
— 12,600 km s^^ on day 7. A "double peaked" emission profile for Hq, H/3, and [Ca II] 
AA7291, 7324 on day 161 is observed, which first appeared between days 95 and 124 and 
persists through at least day 333, the latest epoch obtained thus far. 

There are several mechanisms that can produce SN polarization in addition to asphericity 
of the electron-scattering atmosphere, including scattering by dust (e.g., Wang & Wheeler 1996), 
asymmetrically distributed radioactive material within the SN envelope (e.g., Hoflich 1995), and 
aspherical ionization produced by hard X-rays from the interaction between the SN shock front 
and a nonspherical progenitor wind (Wheeler &; Filippenko 1996). While subtle combinations 
or fine tuning of models using these production mechanisms cannot be ruled out, the primary 
polarization characteristics seen in SN 1999em are exactly what one expects from an aspherical, 
electron-scattering, expanding atmosphere: a wavelength-independent continuum polarization with 
polarization peaks at the locations of P-Cygni troughs and depolarization across strong emission 
features. From the models of H91, we infer that SN 1999em thus had a global asphericity of 
> 7% and a symmetry axis that remained fixed throughout the first 161 days of its photospheric 
development. The observed temporal increase in continuum polarization for SN 1999em may in- 
dicate greater asphericity closer to the explosion center. However, the temporal constancy of its 
EPM-derived distance argues against significant geometric change and suggests that the observed 
polarization increase may be due to the changing density structure of the electron-scattering region. 

The low polarization found for SN 1999em and the lack of significant systematic scatter in 
the Hubble diagram for previous SNe II EPM distances suggest that significant asphericity at 
early times is not prevalent among this SN class. This conclusion is encouraging for cosmological 
applications of SNe II, and is in accord with the observed trend that SN asphericity decreases 
with increasing envelope mass (Wheeler 2000). The anticipated increase in polarization as the 
photosphere recedes through the hydrogen envelope is indeed observed, although the increase is 
quite modest and the interpretation is hampered somewhat by the changing density structure of 
the SN photosphere and envelope. However, imprinting even a small degree of asphericity on the 
outer regions of the thick hydrogen envelope of a SN II-P requires a very large explosion asymmetry, 
since asymmetric explosions tend to turn spherical as they expand (Chevalier & Soker 1989; Hoflich, 
Wheeler, &: Wang 1999). Whether the inferred asymmetry necessitates explosion asymmetry as 
extreme as the "bipolar," jet-induced mechanism proposed by Khokhlov et al. (1999) will have to 
await detailed modeling and acquisition of more spectropolarimetry of core-collapse events. 

This study brings the total number of SNe studied in detail with spectropolarimetry to four. 
Clearly, additional spectropolarimetry of core-collapse SNe, and SNe II-P in particular, are needed 
to fully assess the nature and ultimate cause of SN asphericity at early times. Nevertheless, it is not 
unreasonable to speculate that explosion asymmetry, or asymmetry in the collapsing Chandrasekhar 
core, plays a dominant role in the explanation of pulsar velocities (e.g.. Burrows & Hayes 1996), the 
mixing of radioactive material seen far out into the ejecta of young SNe (e.g., SN 1987A; Sunyaev 
et al. 1987; Lucy 1988; Tueller et al. 1991), and even gamma-ray bursts (MacFadyan & Woosley 
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1999; Wheeler 2000; Wheeler et al. 2000). With many large telescopes with polarimeters currently 
(or soon to be) on line (e.g., Keck, VLT, Gemini), prospects for testing these connections are bright. 
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A. APPENDIX 

A.l. DETAILS OF THE REDUCTIONS 

Miller et al. (1988) discuss the reduction of data obtained using the Lick spectropolarimeter, 
and similar procedures are followed for data acquired with the Keck spectropolarimeter. Here 
we focus on three issues specific to the data acquired for SN 1999em: determining the Stokes 
parameters with only 3 of the usual 4 waveplate position observations, instrumental polarization, 
and sources of systematic uncertainty in SN spectropolarimetry. 

For a given observation, both the Lick and Keck spectropolarimeters place two parallel, per- 
pendicularly polarized beams (B and T for "bottom" and "top") onto the CCD chip, with the 
polarization angle selected by means of a rotatable, achromatic half-wave plate. Rotating the wave 
plate by an angle rotates the plane of polarization by 26. Thus, a 45° rotation exchanges the 
positions of the two parallel beams on the chip. Generally, four observations are required ("1," 
"2," "3," "4"), corresponding to the four waveplate positions (typically observed in this order) 
0°, 45°, 22.5°, and 67.5° (polarization angles 0°, 90°, 45°, and 135° in the plane of the sky), in order 
to measure q and u with sufficient redundancy in the data that changes in seeing, transparency, 
or instrumental response will not affect the results. The main purpose of swapping the location 
on the CCD of the two beams in successive exposures is to remove differences in the instrumental 
response between the B and T beams from the calculation of the Stokes parameters. It is possible 
to compute a gain ratio, gi, which describes the relative response of the B and T beams for the 
first two waveplate positions: 

A gain ratio 52 for the second pair of observations can be similarly calculated. In general, we 
have found gi ~ 52- In order to calculate q and u with the day 159 data for SN 1999em that 
lacked the fourth exposure (waveplate position 67.5°), we had to assume that gi = g2, and use 
gi to remove instrumental response differences between the top and bottom beams in the third 
exposure. While we were thus able to form the complete set of Stokes parameters from this limited 
data set, increased potential for systematic errors obviously existed. Fortunately, we obtained an 
additional (complete) set of observations just four days later at Keck, which largely confirmed the 
earlier results (Fig. 20). (Unfortunately, since the zero point of the spectrograph shifted by ^ 4° 
between the two observations, we could not directly use the fourth waveplate position from day 
163 to reduce the day 159 data without mixing the Stokes parameters; the gratings employed also 
had very different resolutions.) For the analysis of late-time SN 1999em polarimetry presented in 
§ 4.1.3, we thus took the weighted average of these two data sets after first forming the Stokes 
parameters for each separately. 

In general, observations of null polarization standards were null to within 0.1%, and the mea- 
sured polarization of polarized standards agreed well with cataloged values, implying little in- 
strumental polarization over the optical passband at both Keck and Lick. On day 7, however. 
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Fig. 20. — Comparison of the polarization across the Ha profile measured on days 159 and 163. 
For clarity, the lo" statistical error bars are only shown for the day 163 data. The generally good 
agreement, particularly obvious in regions of high signal-to-noise ratio, suggests no large systematic 
error associated with the problematic reduction of the day 159 data. 
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Fig. 21. — Observed polarization measured on 1999 November 5 for SN 1999em and the null 
standard, HD 12021, shown (a) before and (b) after removal of instrumental polarization. 
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there was clear evidence for systematic problems with all of the data collected by the blue arm 
of the spectrograph: both our object and null standards showed a rise of up to ~ 0.4% at blue 
wavelengths (Fig. 21a), and polarization standards showed changes in polarization angle at blue 
wavelengths that disagreed with published values. Martel (1996) found similar effects with the 
Lick spectropolarimeter, and concluded that the most likely source was instrumental polarization 
caused by defects in the aluminum and/or overcoat layers of the primary and secondary mirrors. 
A similar situation has been reported at the Palomar 5-m Hale telescope (Ogle et al. 1999). The 
only null standard observed on day 7 was HD 12021, a star flagged by Schmidt et al. (1992b) as 
possessing some intrinsic polarization, and we did not observe any other null standards with this 
setup until after the primary was realuminized on 2000 February 14. Unfortunately, realuminiza- 
tion alters the polarization characteristics of the mirrors (previous dates for realuminization of the 
Lick 3-m primary are 1996 February 6, 1992 December 8, 1990 February 5, and 1987 March 16; D. 
Severinsen 2000, private communication). Thus, we had to rely on HD 12021 to correct the data 
from this night. 

As in previous studies, we assume that the telescope's intrinsic polarization is produced by the 
mirrors, although polarization due to the dichroic is also possible (see Martel 1996 for a discussion 
favoring the mirrors as the likely cause). To remove instrumental polarization, it is important to 
subtract it as the first step in the polarimetric reduction, before any other corrections are made. 
Since the Lick 3-m is an equatorially mounted telescope with its primary and secondary always 
maintaining the same orientation with respect to absolute direction in the sky, it is necessary to 
carry out the instrumental polarization subtraction with both the instrumental and object polariza- 
tion determined for the same position angle (P.A.) of the spectrograph. That is, for an object with 
observed (i.e., unrotated, "raw") Stokes parameters (9obs7^*obs), the corrected Stokes parameters 
iqcoT,Ucov) are given by 

9cor ^ 9obs ~ ftnst) (A2j 

''^cor ^ ^obs ~ ^inst- 1-^3 j 

Here (^instj ^inst) are the instrumental polarization rotated to the object's P.A. by 

9inst = fcuu cos{2APA) + Unuu sin(2AP^), (A4) 

Uinst = -^nuu sm{2APA) + Unuw cos{2APA) , (A5) 

where (gnuih ^inuu) are the observed, unrotated Stokes parameters for the null standard, and APA = 



Note that for a telescope like Keck-I, which has an altitude-azimuthal mounting, the mirrors maintain their 
orientation relative to the horizon. That is, if instrumental polarization exists at Keck, two separate observations of 
a null standard will produce identical (gnuii, Wnuii) if both are made with the spectrograph at the same orientation 
relative to the horizon, whereas at Lick they will be identical if both are made with the same spectrograph orientation 
relative to absolute direction in the plane of the sky. To correct Keck data one could define APA in equations (A4) 
and (A5) as APA = {[PAohj — f ^opj_par] — [f^nuii — ^^nuii_par]), whcrc "par" refers to the parallactic slit P.A. 
that minimizes differential light loss (i.e., is perpendicular to the horizon; see Filippenko 1982). 
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To derive the instrumental polarization from our observation of HD 12021, we assumed that 
the star's intrinsic polarization was accurately recorded by the red arm of the spectrograph, a belief 
supported by the lack of systematic trends seen in the P.A. of the polarization standards on the red 
side. Therefore, we first removed ~ 0.2% from the measured polarization of HD 12021, presumably 
leaving us with just the "instrumental" portion. We fit low-order splines to the remaining instru- 
mental Stokes parameters rotated to P.A. = 0° (Fig. 22), and then subtracted this spline as the 
first step in the polarimetry reduction from all objects observed on this night (after first rotating 
the object's observed polarization to a P.A. of 0°). After correction, the polarization of SN 1999em 
became flat at blue wavelengths (Fig. 21b), and the polarization angle of the polarization stan- 
dards became more stable, suggesting that we properly removed the instrumental portion of the 
polarization. 

One final issue is the degree systematic errors play in the determination of SN polarization. 
To test the stability of the Lick polarimeter, we took five complete sets of observations in suc- 
cession on 2000 July 28 of the polarization standard HD 19820, without rotating the slit between 
observations, and computed the observed y-band polarization for each set. Although the statis- 
tical uncertainty of each individual broadband measurement was very low (cg^u < 0.008%), we 
found {qv) = 3.122% ± 0.040%, (uy) = -3.663% ± 0.015%, where the quoted uncertainties in 
{qv) and {uy) are the standard deviations of the five individual measurements (note that the re- 
sulting pv = 4.81% ± 0.03% compares well with the published value pv = 4.79% ± 0.03% [Schmidt 
et al. 1992b]). Clearly, systematic uncertainties dominate over statistical ones at the 0.05% level. 
Note, however, that this uncertainty was never seen to manifest itself as artificial line features; 
rather, slight shifts in the overall continuum level were responsible for the differences. 

Barring intrinsic polarization changes in SN1999em over the timescale of a few hours. Figure 23 
demonstrates that systematic errors also dominate statistical uncertainties for SN 1999em. In the 
worst case (day 49, q parameter), successive exposures disagree by more than 9a (Fig. 24). In 
addition to the systematic instrumental uncertainty described earlier, another error source for SN 
spectropolarimetry is the light contributed by the host galaxy, either in the background regions 
or mixed into the SN spectrum. Broadband imaging polarimetry of the diffuse light in nearby 
spiral galaxies indicates polarizations of several percent (e.g., Scarrott, Rolph, & Semple 1990; 
Sholomitskii, Maslov, & Vitrichenko 1999; see also Brindle et al. [1991] in which the central 6" of 
NGC 1637, the host of SN 1999em, is found to be polarized at 1.34% in the B band). While host 
galaxy light is typically only a few percent as bright as the SN, removal of different background 
regions from a SN spectrum in successive sets of exposures (as can happen when the P.A. is rotated 
to maintain the parallactic angle) could certainly affect the overall polarization level at the 0.1% 
level. One possible remedy is to maintain a constant slit P.A. for all exposures, regardless of 
the parallactic P.A. However, this risks allowing polarized galaxy light to be mixed in with the 
SN light itself when observing far away from the parallactic angle, which can produce a spurious 
wavelength dependence to the continuum polarization (e.g., Barth et al. 1999). Variable seeing and 
poor guiding may also play a role. Such considerations may explain the differences seen between 
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Fig. 22. — Instrumental polarization removed from data obtained on 1999 November 5 with the 
Kast double spectrograph at the Lick 3-m telescope. Shown here is the spline fit {thick line) to the 
instrumental q and u Stokes parameters rotated to a telescope P.A. of degrees. The instrumental 
polarization was derived using the null standard HD 12021 on this night; see text for details. 
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Fig. 23. — Evolution of the broadband polarization of SN 1999ein from day 7 to day 161 after 
discovery, displayed as flux- weighted averages over 5050 — 5950 A (approximates observed V band) ; 
due to its different spectral coverage the value for day 7 was computed over 6050 — 6950 A). For 
each epoch, the two individual observations (small open symbols) are connected to the combined 
value (large filled symbols) with a dotted line. Error bars represent the la statistical uncertainties. 
Systematic effects clearly dominate the uncertainty in the broadband polarization measurements. 
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Fig. 24. — Stokes q parameter measured by successive observations of SN 1999em on 1999 December 
17. The second observation appears to be shifted by ~ 0.2% with respect to the first. The offset 
may arise from polarization differences in the host galaxy light in the background regions, due to 
the P.A. rotation made between the two observations. Fortunately, this systematic effect does not 
appear to influence the relative strength of polarization modulations across line features, such as 
the dip seen here at Ha. 
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successive exposures of SN 1999em. Fortunately, while the overall polarization level may be affected 
by diffuse host galaxy contamination, sharp line features are not altered (see Figs. 24 and 20). 
Nevertheless, assigning a la systematic uncertainty to SN continuum polarization measurements 
of ~ 0.1% seems warranted, especially for SNe in dusty spiral galaxies. 



A. 2. Interpretation of Spectropolarimetric Data With Lo^v Polarization 

Studying an object possessing nearly zero continuum polarization but sharp polarization 
changes across line features demands careful attention to detail in the polarimetric analysis. As an 
aid to understanding some of the subtleties involved, we present the following analysis of artificial 
data. 

Figure 25 shows "observed" spectropolarimetry for an object with a polarized continuum 
(^cont = 0%,Ucont = 0.2%), oue polarized line centered at 5700 A {qpoi = 0.6%), Upoi = 0.2%), 
and one unpolarized line centered at 5300 A [qunpoi = O.OVo, Uunpoi = 0.0%), to which we have 
added interstellar polarization (qisp = O-OYo, uisp = —0.5%) with a "Serkowski curve" (Serkowski 
1973; Wilking, Lebofsky, & Rieke 1982) peaking at 5500 A. Gaussian noise oi aq = 0.2%, au = 0.1% 
has also been added to the original q and u data. For ease of analysis, note that all the ISP is 
artificially in the u direction and all the excess line polarization is in the q direction. 

To begin, the unpolarized line is used to estimate the ISP contaminating the object's polariza- 
tion (i.e., § 4.1). This is done by calculating the flux- weighted and continuum-subtracted averages 
of q and u across the unpolarized line's profile, with qisp found by measuring 

_ E(/g) - ifQ)cont] ..^s 

qiine — isr^/r r M ' V-^OJ 

[l^U - Jcontn 

and setting qisp = qune, where {fq)cont and fcont represent fits to the continuum level underlying the 
line in the Stokes parameter flux and total flux, respectively. A similar procedure is carried out to 
determine ujsp. The fit to {fq)cont is preferably taken as the interpolated value of the means found 
for background regions on either side of the line profile. If the background is complicated, however, 
the continuum may be drawn by hand under the line feature in the fq spectrum, although the exact 
placement may be quite subjective. The continuum fits to the Stokes parameter fluxes for the line 
region 5270 — 5330 A are shown in Figure 26; background regions were chosen to be 5150 — 5250 A 
and 5350 — 5450 A. Using equation (A6), we determine qisp = —0.10% ± 0.07% and ujsp = 
—0.51% lb 0.04%, in good agreement with the actual values. To translate this into p and 6 we use 
the "debiased" estimator of p (see below) to obtain pjsp = 0.51% ± 0.04% at 6* = 130° ± 4°, in 
good agreement with the known value. 

The resulting "intrinsic" spectropolarimetry is shown in Figure 27; in order to simplify the 
following analysis by keeping all of the extra line polarization directed parallel to the q axis, we 
have removed the actual, not the measured, ISP. Translating q and u into p can now proceed by 



54 



4.5 
3.0 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

(a) II Example Data 

(with ISP) 



1.5- 



0.0 



^ 



^ 



0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 

0.4 
0.2 
0.0 
0.2 
0.4 
0.6 



s -0.2- 



- (c) 



Unpol. line 

^— M 1 I I 



Pol. line 

I I I 




wrt^ 



\—\—\ — \—\—\ 




5000 



J L 



1 



J I L 



1 



J L 



5200 5400 5600 

Wavelength (A) 



5800 



■friH h-h-F- 




J L 



6000 



Fig. 25. — Artificial data, showing (a) the "observed" flux, {b, c) q, u, and (d) position angle {6) 
in the plane of the sky, uncorrected for the effects of ISP. The "noise-free" Stokes parameters are 
overplotted {thick line) in (b) and (c). 



55 



1.0 
0.5 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Example Data -|- ISP Determination 




5000 5100 5200 5300 5400 

Wavelength (A) 



5500 



5600 
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calculating any of 

(A7) 

(A8) 

(A9) 
(AlO) 

where pdeb assumes the sign of (g^ + u^ — [a'g + cr'i]) and popt is formally defined only when ptrad > 
o'iptrad)- RSP is found by rotating the q and u coordinate system through an angle 9 {6 = 
(1/2) arctan[u/q]), which places all of the polarization along the rotated q parameter. A smoothed 
version of 9 is always used to form the RSP (Fig. 27d). Some of the theory underlying Pdeb-,Popt-, 
and RSP is given by Stockman & Angel (1978; Pdeb)-, Wang, Wheeler, & Hoflich (1997; Popt) and 
Trammel, Dinerstein, & Goodrich (1993; RSP). A discussion regarding the biases introduced by 
using either ptrad or Pdeb on data with low p/ap is presented by Miller et al. (1988). 

Shortcomings inherent to all four estimators of p are apparent in Figure 28. The "traditional" 
polarization estimator, ptrad-, is biased high, and does a poor job detecting the unpolarized line 
feature at 5300 A. The high bias results from its positive-definite definition (eq. [A7]) as well as the 
positive skewness exhibited by its probability distribution (e.g.. Miller et al. 1988), with the former 
cause dominating for the low p/cTp used in this example. Note that Ptrad becomes a progressively 
worse estimator of p as p/cTp decreases, placing a cautionary note on interpretation of increases 
in Ptrad seen in spectral regions where the signal-to- noise ratio is low (e.g., line troughs). The 
"debiased" polarization estimator, Pdeb-, also poorly represents the true polarization (Fig. 28b), since 
its probability distribution is double-peaked and results in a preponderance of negative polarizations 
in many pixels whenp/(Tp is low (Miller et al. 1988). The "optimal" polarization estimator (Fig. 28c) 
has fewer large negative data spikes (formally, popt is undefined at these locations) and overall gives 
a better representation of the true polarization than either ptrad or Pdeb- 

Figure 28d shows that RSP yields an accurate representation of the continuum polarization 
level, with a well-behaved Gaussian distribution of values about the true continuum value. In 
places where 6 changes rapidly (e.g., across the polarized line feature at 5700 A), however, RSP 
underestimates the true polarization since 9 smooth 7^ 6* in the line (Fig. 27d). When this happens, 
some of the line's polarization will be contained in the rotated u parameter (URSP). Thus, when 
determining line polarization with RSP, it is often necessary to inspect both RSP and URSP. 
Figure 29 shows the data for both the unpolarized and polarized lines, binned by 10 A to reduce 
noise. Since 9 does not change across the unpolarized line, RSP accurately estimates PunpoU showing 
only random residuals about in URSP. Conversely, RSP underestimates Ppoi , with consecutive bins 
consistently showing large deviations (all with the same sign) from zero in URSP. It is generally not 
prudent to rectify this situation by following the 9 curve more closely (i.e., smoothing it less), since 
RSP then becomes subject to the same biases that affect Ptrad- Indeed, in the limit 9smooth -^ d, 
RSP -^ ^/^T^. 
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Fig. 27. — Artificial data with ISP removed. The "true" values of the Stokes parameters are 
overplotted (heavy line) in (6) and (c), along with the smoothed fit (smoothed and median filtered 
with a boxcar of width 100 A, our typical smoothing value) to the polarization angle (6) used to 
determine the RSP, in (d). Note that the sharp polarization angle change across the polarized line 
at 5700 A is substantially missed by the fit. 
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Fig. 28. — Four estimates of the polarization of the artificial data, with the "true" (noise-free) 
polarization overplotted {thick line), (a) The "traditional" p is biased high; (b) the "debiased" 
p shows many large negative values; (c) the "optimal" p is better behaved, but still shows a few 
large, negative spikes; (d) the rotated Stokes parameter provides a good representation of the 
true continuum polarization, although the polarization of the line feature at 5700 A is somewhat 
underestimated (see Fig. 29). 
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Fig. 29. — (o) Normalized flux, (&) RSP, and (c) URSP for the artiflcial data in the vicinity of the 
unpolarized hne feature at 5300 A and the polarized line feature at 5700 A. The true polarization 
level is overplotted (thick line) in (b). The fact that the RSP underestimates the polarization in 
the polarized line could be deduced by inspection of the line's URSP, which shows large residuals 
with a consistent sign in the three central pixels. 
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There is no perfect estimator of the true polarization for data with low p/cTp, and the best 
one to use generally depends on the polarization properties of the type of data being analyzed. 
For data with p/(Jp ^ 1 (e.g., SN 1998S), all four estimators produce essentially identical results. 
When p/cTp < 3, RSP is most robust, but only when 9 varies slowly with wavelength: sharp 9 
changes necessitate examination of URSP as well. Alternatively, one can use a coarser binning 
to make p/cTp > 3 so that pdeb and Popt become reliable, but care must then be taken to ensure 
that polarization features are not altered by the decreased resolution (i.e., the binning should 
be significantly narrower than the scale over which polarization features appear). Finally, when 
summing over a very large wavelength band to determine an average polarization (e.g., pv)-, taking 
the flux-weighted average of unbinned data (or data with fine binning) and then computing either 
Pdeb or Popt is the method of choice (see, e.g.. Miller et al. 1988). In this paper, we use pdeh to 
measure broadband polarization and RSP to characterize the spectropolarimetry. 

For the artificial data, we can now estimate the line's polarization from the ISP-subtracted 
data using equation (A6) and fitting the Stokes parameter fluxes as before. For the polarized line, 
we measure qpoi = 0.58% ± 0.07% and Upoi = 0.21% it 0.04%, calculated over the line region 
5670 — 5730 A with background regions 5550 — 5650 A and 5750 — 5850 A. For the unpolarized 
line we measure Qunpoi = —0.01% ± 0.08% and Uunpoi = —0.01% ± 0.04% using the same line 
and background regions described earlier for determining the ISP. Our analysis techniques have 
successfully recovered the input values, though the uncertainty in the ISP removal would add 
systematic uncertainty to the error budget. 
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